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We have shown within the quasistatic approximation that the giant fluctuations of local electro- 
magnetic field in random fractal aggregates of silver nanospheres are strongly correlated with a local 
anisotropy factor S which is defined in this paper. The latter is a purely geometrical parameter 
which characterizes the deviation of local environment of a given nanosphere in an aggregate from 
spherical symmetry. Therefore, it is possible to predict the sites with anomalously large local fields 
in an aggregate without explicitly solving the electromagnetic problem. We have also demonstrated 
that the average (over nanospheres) value of S does not depend noticeably on the fractal dimension 
D, except when D approaches the trivial limit D — 3. In this case, as one can expect, the average 
local environment becomes spherically symmetrical and S approaches zero. This corresponds to the 
well-known fact that in trivial aggregates fluctuations of local electromagnetic fields are much weaker 
than in fractal aggregates. Thus, we find that, within the quasistatics, the large-scale geometry does 
not have a significant impact on local electromagnetic responses in nanoaggregates in a wide range 
of fractal dimensions. However, this prediction is expected to be not correct in aggregates which 
are sufficiently large for the intermediate- and radiation-zone interaction of individual nanospheres 
to become important. 



I. INTRODUCTION 

Electromagnetic properties of fractal nanostructures 
have continuously attracted attention since the late 1980- 
ies due to their rather unusual physical properties and 
the possibility of numerous applications, as described 
in several reviews of the subject Q, S S S S 1 . Of 
special interest are aggregates of metal nanoparticles in 
hydrosols and percolation clusters (metal-dielectric com- 
posites) which have, in particular, exhibited the effects of 
giant enhancement of nonlinear-optical responses 0,0,0, 
[l0| , inhomogeneous localization of electromagnetic eigen- 
modes [Tllll^ . and optical memory [lH HI 111 HI fl7| . 

Theoretical description of the electromagnetic re- 
sponses of disordered fractal aggregates has been closely 
interconnected with numerical simulations. This is due 
to the fact that a fully analytic solution to the problem of 
interaction of an electromagnetic field with a large ran- 
dom fractal aggregate has not been devised. Some ap- 
proximate theoretical approaches were based on the first 
Born [Tsj and mean- field |l9j| approximations, approxi- 
mations based on few-body interaction (binary [2(1 l2lj 
or binary-ternary 01 approximations) , and various phe- 
nomenological scaling laws H |U H3 El |H H . The 
first Born and the mean-field approximations are not, 
generally, useful in the spectral regions where excitation 
is resonant. While off- resonant electromagnetic proper- 
ties of fractal a ggr egates are of interest in atmospheric 
physics [2{J l26l l27j , the research in electromagnetics of 
metal fractal aggregates is primarily focused on resonant 



interactions. The few-body approximations and the scal- 
ing laws proved to be very useful for qualitative theoret- 
ical description at the early stages of research. However, 
increasingly more realistic simulations revealed that these 
approaches do not provide quantitative results. Cur- 
rently, they are effectively obsolete. A brief overview 
of the progression of numerical models used to simulate 
electromagnetic responses of fractal aggregates is given 
in the next paragraph. 

The theoretical and computational description has 
been primarily based on a model of an aggregate of 
N touching identical spherical nanoparticles. Each 
nanoparticle, and sometimes the aggregate as a whole, 
are assumed to be much smaller in size than the ex- 
ternal wavelength. (Polydisperse aggregates built from 
spheres of different size have also been recently ad- 
dressed [^IHIjicj.) In order for an aggregate to be con- 
sidered fractal, the number of primary spheres, must be 
large, typically, > 10 3 . Initially, simulations were based 
on the dipolc approximation. In this approximation, each 
sphere is assigned a dipole moment located at its center. 
The spheres then interact with each other and the exter- 
nal field via dipole radiation fields as described by 37V 
coupled-dipole equations |2l|. In the late 1980-ies and 
early 1990-ies, numerical solution of dense linear systems 
of > 10 3 equations was a difficult computational task. 
Therefore a model of diluted aggregates was adopted and 
used, for example, in Refs. g Ell El EJ |H H|. 
According to this model, an aggregate of N touching 
spheres (where N can be very large) is diluted, i.e., 
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spheres were randomly removed from the aggregate with 
the probability 1 — p, where p <C 1. Then the coor- 
dinates of the remaining spheres arc rescalcd according 
to r — > p x l D v, where D is the fractal dimension. This 
procedure does not change the density-density correla- 
tion function of the aggregate in some intermediate re- 
gion. However, it does change the local structure of 
the aggregate substantially. The few-body approxima- 
tions and scaling laws were largely validated with the 
model of diluted aggregates. However, when computa- 
tions with non-diluted clusters became feasible, it was 
found that both the few-body approximations and the 
scaling laws are inaccurate |35j . The deviation from the 
scaling laws has been explained by the phenomenon of 
inhomogeneous localization however, the theoretical 
relation of this phenomena to the aggregate geometry has 
not been clarified. Additionally, it has been well known 
that account of excitation of higher multipole modes is 
important for touching nanoparticlcs, even when the size 
of each nanoparticle is much smaller than the external 
wavelength [3g, |33, l2a> l^M ■ I n particular, the dipole ap- 
proximation failed to properly describe experimentally 
observed red shifts in extinction spectra of colloid aggre- 
gates [3^, 0] . To remediate this problem, a phenomeno- 
logical model of geometrical renormalization have been 
introduced [35l l4lj and, recently, computations beyond 
the dipole approximation have been performed |42| . 



The combination of findings contained in the above- 
cited references strongly suggest that the local structure 
of aggregates is of primary importance. However, the lo- 
cal structure of random fractal nanoaggregates has not 
been so far the focus of research. In this paper we, for 
the first time, point to a strong correlation between the 
anisotropy of local environment and enhancement of lo- 
cal field in fractal aggregates within the quasistatic ap- 
proximation. In particular, we find that the correlation 
coefficient of the local anisotropy factor S (introduced be- 
low) and the value of a local squared dipole moment |d| 2 
can be as high as 0.75 and tends to grow with the wave- 
length. We have found that the average local anisotropy 
factor is almost independent of fractal dimension in the 
range 1.7 < D < 2.8. Note that this result is expected 
to change in large aggregates where intermediate- and 
far-zone interaction is important. 



The paper is organized as follows. In Section UTI the 
local anisotropy factor is introduced. The dependence of 
the local anisotropy factor on the fractal dimension of 
aggregates and other parameters for computer-simulated 
fractals is discussed in Section lTTTl Section llVl contains re- 
sults concerning the correlation of local electromagnetic 
fields and the local anisotropy factor. The electromag- 
netic calculations in this section were performed with the 
method of coupled multipoles [33,|42|, e.g., without the 
dipole approximation. Finally, Section[V]contains a sum- 
mary of obtained results. 



II. DEFINITION OF THE LOCAL 
ANISOTROPY FACTOR 

The definition of local anisotropy factor introduced in 
this paper is based on an analogy with ellipsoids. An el- 
lipsoid is a geometrical object that can exhibit cither per- 
fect spherical symmetry, or strong anisotropy, depending 
on its eccentricity. 

Consider a general ellipsoid excited by a linearly polar- 
ized monochromatic external wave of amplitude Eo. In 
the quasistatic limit, the polarization P inside the ellip- 
soid is independent of position and can be found from 
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where the tensor Q is given by 



En 



(1) 



Q 



Go(0,r')dV 



(2) 



Here Go(r,r') is the regular part of the quasistatic free- 
space dyadic Green's function for the electric field. The 
integral is taken over the volume of the ellipsoid, V, and 
is independent of position. Therefore, it is evaluated at 
the center of ellipsoid, r = 0. A unique property of ellip- 
soids is that Q is diagonal in the reference frame whose 
axes are collinear to the main axes of the ellipsoid. Corre- 
spondingly, if E a are the Cartesian components of the 
external electric field in the same reference frame, the 
solution to is 
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where v a are the depolarization factors related to the 
principal values of Q by 



1/3 - Q c 



(4) 



In the case of spherical symmetry (e = 0), Q a = and 
v a = 1/3. For an ellipsoid of nonzero eccentricity, the 
depolarization factors become different from 1/3. Thus, 
for example, if e = 1, we have v\ = vi = 0, ^3 = 1 
for an oblate ellipsoid (infinitely thin circular disk) or 
V\ = V2 = 1/2, ^3 = for a prolate ellipsoid (infinitely 
thin needle). The anisotropy factor S can be defined as 
dispersion of the depolarization factors: 



S 2 = 



(5) 



Obviously, this parameter is zero for a sphere and positive 
for any ellipsoid of nonzero eccentricity. In particular, for 
the infinitely thin needle, S = l/3\/2 and for an infinitely 
thin circular disk, S = V2/3. The latter is the maximum 
possible value for S given the constraint J2 a v a = ^- 
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Now we extend the definition of the depolarization ten- 
sor to include particles of arbitrary shape. Namely, for 
an arbitrary system occupying some volume V, wc define 



z>(r) 



G (r,r')dV 



(6) 



where / is the unity tensor. If V is of general shape, 
the result of integration in the right-hand side of ||BJ| is 
position-dependent. Therefore, the tensor u(r) depends 
on the point r it is evaluated at and is referred to here as 
local. Similarly to the case of ellipsoids, this tensor can 
be diagonalized. Then we can use the principal values 
u a (r) to calculate the anisotropy factor according to JSJ. 

In this paper we consider aggregates of (possibly, poly- 
disperse) spheres whose centers are located at points rj 
and radii are denoted by a;. In this case, the expression 
for i>(r) is simplified. We use 
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G (r,r')dV 
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to obtain 



«iG (r,ri) , if |r - n\ > a, t 
, if jr- rij < dj 

(7) 



Vi = z>(r 4 ) = \l~^2 VjGo(ri, rj) 



(8) 



Here = 47raf/3 is the volume of the i-th sphere and 
the components of Go(i"j, r j) are given by 



principal values. Second, the principal values v a (r) are 
not constrained, in general, by the conditions < v a < 1 
and v a = 1. This also distinguishes them from the 
Bergman-Milton spectral parameters. Next, the param- 
eter Si depends on the coordinates of all nanoparticles in 
the aggregate with j ^ i. However, due to the fast cubic 
decay of the near-field component of the dipole radiation 
field, the neighbors within few coordinate spheres of the 
i-th site give the largest input to Si. This justifies the lo- 
cality of Si, as it only weakly depends on the large scale 
structure. This statement needs to be qualified in ag- 
gregates large enough so that interaction in the far-zone 
becomes important. Even without account of retarda- 
tion, the locality of Si can be violated in aggregates with 
the fractal dimension close to 3 (or in random non-fractal 
composites), due the logarithmic divergence of the inte- 
gral / r~ 3 d 3 r at infinity. We do not expect these effects 
to be important in most aggregates of practical interest 
with the fractal dimension in the range D < 2.7 and do 
not consider them in this paper. 

Finally, the introduced parameter is not sensitive to 
the wavelength and electromagnetic properties of the 
scattering material. Therefore, we do not expect it to be 
a good indicator of local electromagnetic response at all 
wavelengths. It is also independent of the incident polar- 
ization. A possible definition of a polarization-sensitive 
anisotropy factor is 



S;(E ) 
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Another possible definition is 
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Sf (E„) = 
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where = r,; — rj and = 

Diagonalization of the tensor i>i and calculation of 
the dispersion of its principal values gives the local 
anisotropy factor Si- This parameter quantifies the de- 
gree of anisotropy of the local environment of the i-th 
sphere. 

A few notes about the introduced definition must be 
made. First, the principal value v a obtained as described 
above are purely geometrical characteristics of an object. 
They are related to the Bergman-Milton spectral param- 
eters 01 only in the special case of ellipsoidal (more gen- 
erally, spheroidal) shape of V. Obtaining the Bergman- 
Milton spectral parameters requires diagonalization of 
the integral operator W with the kernel Go(r, r'), r, r' e 
V. This is a much more complicated task than diagonal- 
ization of the tensor Q(r) = J v Go(r, r')d 3 r' at a given 

point r. In particular, Q(r) is three-dimensional, while 
W is infinite-dimensional. Correspondingly, the number 
of Bergman-Milton parameters is infinite (although only 
three of them have non-zero oscillator strengths in the 
case of spheroids), while the tensor v{r) has only three 



Note that the definitions l|10 [l .l)ll [l are not used in this 
paper. 



III. RESULTS: SIMULATIONS OF 
GEOMETRICAL PROPERTIES 

Since the unique electromagnetic properties of colloid 
aggregates are often attributed to their fractal struc- 
ture, we have studied computer-generated aggregates 
with various fractal dimensions. We have generated 
quasi-random off-lattice aggregates with varying fractal 
dimension D using the algorithm described in [4^ |. This 
algorithm simulates the stochastic dynamics of individual 
nanoparticles and sub-aggregates in a solution with the 
account of random (Brownian) forces, as well as deter- 
ministic intcrparticle (the Van-dcr-Waals and Coulomb) 
and external potentials. Discrete Newtonian mechanics 
was implemented with a sufficiently small time step, such 
that the spatial translation of any particle (sphere) at 
each step is much smaller than its diameter. Rotation of 
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aggregates was taken into account. We have used both 
monodisperse (<Xi = const) and polydisperse nanosphcrcs 
{a,i were randomly distributed according to the Poisson 
distribution). The fractal dimension of obtained aggre- 
gates was tuned in the interval 1.7<Z3<3.0by varying 
the initial density of spheres prior to the aggregation pro- 
cess. The numerical value of D was calculated from the 
linear regression of the pair density-density correlation 
function which, in the intermediate asymptote region, 
has the scaling form g(r) cx r D ~ 3 . 

The aggregation was simulated in a cubic volume 
with elastically reflecting boundaries. In the limit of 
low initial concentration of particles and the size of the 
cubic cell of ~ 200a or more, the obtained aggregates 
have the typical fractal dimension D w 1.7. When the 
initial concentration increases, D approaches the trivial 
limit D = 3. As a graphical illustration of generated 
fractals, we show in Fig. 1 a large aggregate and val- 
ues of the local anisotropy factor S at some selected sites. 




S=2.03 



S=1.48 



S=0.01 



average S for two types of lattice aggregates traditionally 
used in electrodynamic calculations. The results are 
shown by centered symbols in Fig. 2. 




Fig. 1. Local anisotropy factor S for selected sites in a large 
aggregate with N = 5000 and D « 1.8. 



Fig. 2. Average local anisotropy factor S vs fractal di- 
mension D for 121 random aggregates with N = 800 in 
each. Numerical value of D was computed separately for 
each aggregate from linear regression of the density-density 
correlation function in the intermediate asymptote region. 
The open circle shows the average values of S for lattice 
Meakin aggregates |3 (D ~ 1.8) and the open square shows 
the same value for a set of Witten-Sander aggregates |45l ) 
(D w 2.5). 



Real colloid aggregates are strongly polydisperse. Typ- 
ically, they contain particles of sizes ranging from 5nm to 
30nm plUg]. We have investigated the dependence S(D) 
for several ensembles of polydisperse aggregates with dif- 
ferent ratios of the maximum and minimum sphere ra- 
diuses, a max and a m ; n . We have used a discrete Pois- 
son distribution of particle sizes with the number of 
samples equal to 11 [the discrete step in particle size 
was Aa = (a max — a m i n )/10]. The dependence of local 
anisotropy factor on the fractal dimensions D is shown 
in Fig. 3. Note that no signifact effect doe to the polidis- 
persity was found. 



We start with a discussion of results for monodis- 
perse aggregates, i.e., for aggregates built of identical 
spheres. In Fig. 2, we illustrate the dependence of 
the average (over individual particles in an aggregate) 
value of the S on fractal dimension D. Aggregates with 
1.70 < D < 2.25 are characterized by moderate average 
values of local anisotropy factor, almost independently 
of D. We can argue that such aggregates differ only 
on large scale but have similar local structure. In 
other words, the local environment of each particle 
is, on average, the same, independently of D. As D 
approaches the critical value D = 3, the local anisotropy 
factor quickly drops. This corresponds to the fact that 
trivial (non-fractal) aggregates are characterized by 
almost isotropic local environment and relatively weak 
fluctuations of density. We have also calculated the 



It is interesting to note that the average local 
anisotropy factor does not depend on the distance of a 
given site from the center of mass of the aggregate. This 
is illustrated in Fig. 4. Here we plot the value of S aver- 
aged over all nanospheres within a spherical shell drawn 
arond the center of mass of the aggregate as a function 
of the shell radius (see figure captions for more detail). 

In Fig. 5, we also plot the fraction of particles in an 
aggregate with local anisotropy factor exceeding 60% of 
the maximum value for that aggregate as a function of 
fractal dimension. It can be seen that in typical aggre- 
gates with fractal dimensions of practical interest, only 
small fraction of particles is placed in highly anisotropic 
environment. In Fig. 6, an example of an aggregate is 
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Fig. 3. Local anisotropy factor S factor vs fractal dimensions 
D for polydisperse aggregates with N = 800 and Poisson 
particle size distribution. The ratio of the maximum and 
minimum particle radii is a max /a m i n = 2 (153 random 
aggregates) for curve (a) and dmax/aWn = 3 (297 random 
aggregates) for curve (b). Numerical value of D was com- 
puted separately for each aggregate from linear regression of 
the density-density correlation function in the intermediate 
asymptote region. 
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Fig. 4. Local anisotropy factor S vs the relative distance to 
the center of mass of an aggregate, r/a. The hystogram is 
built with the step _R 9 /10, where R g is the gyration radius of 
the aggregate, and S was averaged over all particles located 
within 10 spherical shells drawn around the aggregate's 
center of mass for N = 10, 000 (a) and N = 3, 000 (b). 
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Fig. 5. Dependence of the fraction AN/N of sites in an 
monodisperse aggregate with the value of local anisotropy 
factor exceeding 60% of its maximum value for the same 
te; N = 800. 




Fig. 6. Sites in a fractal aggregate (D « 1.8, N = 800) with 
relative values of the local anisotropy factor exceeding 80% 
of the maximum value for the same aggregate, S max = 2.29. 




(b) 



(c) 




200 400 600 800 

Nanoparticle index 

Fig. 7. Local anisotropy factor S for different particles in a 
fractal aggregate with iV = 800 and D m 1.8 (a) compared 
to those in random gas of identical hard spheres. Average 
distance between centers of two nearest neighbor spheres 
Rnn/a = 2.12 (b) and i? nn /a = 2.90 (c). 
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shown with the sites of relatively high local anisotropy 
emphasized by different color (shade of grey). 

Finally, we compare the local anisotropy factors for 
all particles of monodisperse fractal aggregate and 
non-fractal random gas of hard spheres (TV = 800 in 
both cases). In Fig. 7a we plot these quantities for a 
fractal aggregate with D as 1.8. All local anisotropy 
factors Si are shown for i = 1, ...,800. In Fig. 7b, c, 
the same quantities are plotted for a random gas of 
identical hard spheres of radius a distributed in a 
volume with the density corresponding to the average 
distance between the centers of nearest neighbor spheres 
equal to R nn ] value of the ratio R nn /ci are indicated 
in the figure caption. It can be seen that the fractal 
aggregate contains sites with much higher values of local 
anisotropy factor than random gas. As one could expect, 
the local anisotropy factors become smaller when the 
density of random gas decreases. However, a fractal 
aggregate, although has zero asymptotic density in the 
limit N — > oo, always retains approximately constant 
fraction of sites with relatively high local anisotropy. 
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Fig. 8. Local anisotropy factor (Si) (thick line) and local 
dipole moments squared |cL| 2 (thin line) for different parti- 
cles in a monodisperse aggregate with N = 150 and fractal 
dimension D — 1.70, D = 1.79, and D = 1.85, computed at 
the wavelength A = 703nm. 
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Fig. 9. Correlation between Si and |cL| 2 as a function of 
wavelength for monodisperse aggregates with TV = 150 and 
fractal dimension D = 1.70, D = 1.79, and D = 1.85. 



IV. COMPARISON OF STRUCTURAL AND 
ELECTRODYNAMIC PROPERTIES OF 
FRACTAL NANOAGGREGATES 

The main idea of this paper is that there is certain 
correlation between local structure and local electromag- 
netic fields in fractal nanoaggregatcs. This assumption 
is confirmed by numerical simulations presented in this 
Section. 

We have computed optical responses of aggregates of 
nanospheres using the method of coupled multipoles |39l 
Calculations were performed for monodisperse ag- 
gregates built of N = 150 silver nanospheres of constant 
radius a = 5nm and placed in vacuum. To facilitate con- 
vergence with the maximum order of multipoles included, 
we have introduced a surface layer of thickness h = 0.05a. 
The dielectric constant of the layer was chosen to be the 
same as that of the vacuum, e = 1. We have used exper- 
imental values of the optical constants of silver with 
finite-size corrections according to |35j . The maximum 
order of the VSHs utilized in the results shown below 
was L = 8. The convergence was verified by control cal- 
culations with L = 16. We note that much larger values 
of L are required for nanospheres in exact contact (h = 0) 
and that the number of the couplcd-multipolc equations 
(with complex coefficients) which must be solved to com- 
pute the optical responses is equal to NL(L + 2). 

In Fig. 8 we plot the quantities Si and (|dj| 2 )/(a 3 i?o) 2 
for three aggregates with fractal dimensions D as 1.70, 
D as 1.79 and D as 1.85, computed at A = 703nm. 
Here (|di| 2 ) is the square of the dipole moment of i-th 
nanosphere averaged over three orthogonal polarizations 
of the external field. Visual correlation of the two curves 
is quite apparent. For a more quantitative estimate, we 
have computed the correlation coefficient r c (S, (|d| 2 )). 
The dependence of r c on the wavelength is shown in Fig. 9 
for three values of fractal dimension and different wave- 
lengths. The maximum degree of correlation is achieved 
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for A = 703nm (0.69 < r c < 0.76). The value of r c de- 
creases monotonously for smaller wavelengths and is in 
the interval 0.49 < r c < 0.54 when A = 505nm. Note 
that the correlation coefficient is expected to increase to- 
wards unity in the spectral region A > 700nm. Also, 
even stronger correlation is expected is a polarization- 
dependent definition of the local anisotropy factor is 
used, such as <|10[) or 1)11(1. Validating these hypothesis 
will be the subject of future work. 

V. SUMMARY 

In this paper we have investigated the statistical cor- 
relation between the local geometrical structure and lo- 
cal electromagnetic responses in fractal aggregates of 
nanoparticles. We have used a realistic aggregation 
model which allows computer generation of quasi ran- 
dom aggregates of variable fractal dimension in the in- 
terval 1.7 < D < 3.0. Electromagnetic calculations were 
carried out using the method of coupled multipolcs, i.e., 
beyond the dipole approximation. 

We have found that the local anisotropy factor S intro- 
duced in Section l^is strongly correlated with the local 
electromagnetic response. For aggregates built of high- 
quality plasmonic materials, the degree of such correla- 
tion tends to increase with the wavelength. The correla- 
tion coefficient between the squared dipole moment of a 
given nanoparticle in an aggregate and a purely geometri- 
cal parameter (local anisotropy factor) reaches the value 
of « 0.75 for A = 700nm. We expect that this correlation 
can become even larger if a properly-defined polarization- 
dependent local anisotropy factor is used and at larger 
wavelengths. 

The introduced parameter S is a universal geometri- 



cal characteristic which can be used for analyzing vari- 
ous complicated aggregates and composites without ex- 
plicit solution of the electromagnetic problem. The dis- 
covered strong correlation suggests that, at least in ag- 
gregates which are small compared to the wavelength, 
the large-scale geometry does not play a significant role. 
Note that in the IR spectral region, subwavelength ag- 
gregates can still be built of hundreds or even thou- 
sands of nanospheres. The IR spectral region is of spe- 
cial interest because of the very low Ohmic losses in sil- 
ver and other noble metals. Correspondingly, heteroge- 
neous nanostructures arc known to exhibit optical res- 
onances of very high quality. This, in turn, results in 
giant amplification of local optical responses. The lat- 
ter phenomenon is currently being actively researched in 
ordered nanostructures, including self-similar chains of 
nanospheres (nanolenses) |4Sl | and long chains of simi- 
lar nanoparticles |49j . Rigorous numerical simulations 
in random nanoaggrcgates are still difficult due to the 
high computational complexity of the associated electro- 
magnetic problem. The introduced parameter S and the 
discovered correlation of this parameter with local elec- 
tromagnetic field allows one to make qualitative predic- 
tions about the sites where the electromagnetic energy is 
localized by very simple means, e.g., without solving the 
electromagnetic problem. 
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